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We present a method to solve the problem of Rashba spinorbit coupling in semiconductor quantum 
dots, within the context of quasi-exactly solvable spectral problems. We show that the problem 
possesses a hidden osp(2, 2) superalgebra. We constructed a general matrix whose determinant 
provide exact eigenvalues. Analogous mathematical structures between the Rashba and some of the 
other spin-boson physical systems are notified. 

PACS numbers: 

INTRODUCTION 

The optical and electrical properties of confined electrons in semiconductor quantum wells, quantum dots and quan- 
tum wires depend on the Rashba spin orbit couplingJlUdJa . The analysis of the behavior of spins in semiconductors 
leads to the construction of new spintronic devices • ^ i s practical to contemplate semiconductor devices 
based on electron spins, because spin is not coupled to electromagnetic noise and hence should have much longer 
coherence times than charge|8|. The spin-orbit interaction in a confined geometry has also interesting consequences 
for the electron spectrum 0]. 

Spin-orbit interactions can arise in quantum dots by various mechanisms related to the electron confinement and 
symmetry breaking. In semiconductors, spin-orbit coupling from the relativistic effect is caused by the electric field 
due to the lack of the inversion symmetry of the certain alloys. Depending on the particular origin of the electric field, 
the spin orbit interaction presents two distinct contributions; these are Dresselhaus and Rashba terms. Dresselhaus 
term is arised from the electric field produced by the bulk inversion asymmetry of the material and Rashba term is 
generated due to the structural asymmetry of the heterostructure. The Rashba splitting has been observed in many 
experiments and it constitutes the basis of the proposed electronic nonostructures. 

In literature, most of the studies about the solution of the spin-orbit effects in quantum dots are carried out 
by means of perturbation theory or numerical methods, because the implementation of the algebraic techniques to 
solve those problems is not very efficient and most of the other analytical techniques do not yield simple analytical 
expressions. The aim of this paper is to provide exact solution of a quantum dot Hamiltonian including Rashba 
and Zeeman term by an algebraic formulation. Analytical solution of the problem has recently been treated by 
employing various techniques 0, 0, 0, flil fl3j |. We demonstrate that the exact eigenf unctions and eigenvalues 
are available for the Hamiltonian of a quantum dot including Rashba coupling in the framework of quasi-exact 
solvability (QES) [3 E3 • As a related topic the concept of quasi exactly solvable systems discovered in 1980s 
has received much attention in recent years, both from the view point of physical applications and their mathematical 
beauty. It turns out that in quantum mechanics there exist such systems; a part of their spectrum can be computed 
using algebraic methods. We also show the corresponding Hamiltonian possesses a hidden osp(2, 2) algebraic structure. 

Our formulation of the Hamiltonian associated with the quantum dots leads to an interes ting consequence: Hamil- 
tonians of the E®e Jahn- Teller 18, 19, 2JJ, |2l|, |22| , two mode bosonic .Tavnes-Cumrmngs |23l 124. l25j and quantum dot 
Hamiltonians with spin-orbit interaction are constructed for similar spin-boson models, though not identical, and can 
be solved with the same mathematical techniques. The connection between the Hamiltonians plays crucial roles to the 
solution and analysis of the Rashba Hamiltonian, because the solutions of the Jahn- Teller and the Jaynes-Cummings 
problems have been focus of interest both from the mathematical and physical point of view, over years and their 
exact solution has been treated by various authors. 

The paper is organized as follows. In section 2, we briefly review the construction of the Hamiltonian that include 
Rashba spin-orbit coupling term. Section 3 is devoted for algebrization of the corresponding Hamiltonian. In this 
section we also discuss the symmetry properties of the Hamiltonian. In section 4, we present a transformation 
procedure that is appropriate to determine the QES of the Hamiltonian. Finally, we conclude our results in section 4. 
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THE MODEL HAMILTONIAN 



The origin of the Rashba spin-orbit coupling in quantum dots is due to the lack of inversion symmetry which causes 
a local electric field perpendicular to the plane of heterostructure. The Hamiltonian representing the Rashba spin 
orbit coupling for an electron in a quantum dot can be expressed as 

Hr = ~J~ (Py ax ~ "PxOy) (1) 

where Xr represents the strength of the spin orbit coupling, and it can be adjusted by changing the asymmetry of the 
quantum well via external electric field and the matrices a x , and a v are Pauli matrices. The Hamiltonian dominated 
by the bulk inversion symmetry term, is known as Dresselhaus spin orbit coupling and it is given by 



Hd = (PyVy ~ Px&x) (2) 



where Ad is the Dresselhaus parameter that depends on the material properties, device design and external electric 
field. In a centrosymmctric crystal like materials it becomes zero, because of the non existence of bulk inversion 
asymmetry. In the case of both Rashba and Dresselhaus interactions exist in the system the Hamiltonian can not 
exactly be solved. During our treatment we have included Rashba term rather than the Dresselhaus term, because 
Rashba term may be dominant, since the typical ratio of the coupling strengths Xr/Xd = 6|ll|. We mention here 
that, the results will only need a trivial modification, when a solo Dresselhaus term is present because Rashba and 
Dresselhaus terms transform into each other under spin rotation. To this end we assume that the electron is confined 
in a parabolic potential 

F=im^ 2 (x 2 +2/ 2 ) (3) 

here m* is the effective mass of the electron and ujo is the confining potential frequency. The Hamiltonian describing 
an electron in two-dimensional quantum dot takes the form 

H = 2ib (P * + P ^ + l 9 ^ + V + (4) 

The term ^gfj,Ba z introduces the Zeeman splitting between the (+)x— polarized spin up and (— )ac— polarized spin 
down. The factors g is the gyromagnetic ratio and \x is the Bohr magneton. The kinetic momentum P = p + eA is 
expressed with canonical momentum p = —ih(d x , d y , 0) and the vector potential A can be related with the magnetic 
field B = Vx A. The choice of symmetric gauge vector potential A = B/2(—y, x, 0), leads to the following Hamiltonian 

1 ( d d \ 1 

-zhu c [y- - x - + -g^Ba z + H R (5) 



where uj c = eB/m* stands for the cyclotron frequency of the electron, ui = y uj 2 + {^f) 2 is the effective frequency. 
From now on we restrict ourselves to the solution of (0) ■ It will be shown that the Hamiltonian © without Dresselhaus 
term is one of the recently discovered quasi-exactly solvable operator 0|. It is well known that the underlying 
idea behind the quasi exact solvability is the existence of a hidden algebraic structure. In the following section we 
obtain an algebraic expression for the Hamiltonian H and discuss its quasi-exact solvability. 



QES OF THE HAMILTONIAN 



The general procedure to solve a differential equation quasi-exactly is to express it in terms of the given Lie algebra 
having a finite dimensional invariant subspace and use algebraic operations. One way to relate the Hamiltonian H 
with an appropriate Lie algebra is to construct its bosonic and fermionic representation. For this purpose let us 
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introduce the following bosonic operators: 



Ah V 4m*w 



(x - iy) + \ j——(d x - id y ) 



H 



° 2 = v~^ (x + iy) + Vi^ a + ^ ) - 

They satisfy the usual commutation relations. The Hamiltonian _ff describing a two-level fermionic subsystem coupled 
to two boson modes can be expressed as: 

H = huj(a^ai + a\a 2 + 1) H ^— (a^ai — a^ a 2) — 











°v 




-1 0\ 


! 


1) 




! 


) 


1 0Q = ( 


v 1 J 



— Ah [(<4 - aiW + (a 2 - af)a^] + -g(iBa (7) 
The Pauli matrices are given by 

/ n 1 \ / n n \ / _1 n \ 

(8) 

Before we turn to explore the solvability of the Hamiltonian H we briefly discuss its Lie (super)algebraic properties. 

The natural step to relate a Hamiltonian with a Lie (super)algebra is to express the Hamiltonian with the generators 
of the relevant symmetry group. One of the major symmetry group candidates for spin one-half particles is the 
supergroup osp(2, 2) which has four even and four odd generators. Its even generators can be represented by bosons 
while odd generators are represented by combinations of the fermions and bosons [2^, [^J. The Lie superalgebra 
osp{2 1 2) can be constructed by extending su(l, 1) algebra whose generators are given by 

J+ = afa 2 , J- = a 2 ai, J Q = - { a f a i + a^ a 2 + l) ■ (9) 
These are the Schwinger representation of su(l, 1) algebra and its number operator is given by, 

N = a^ai - a^a 2 (10) 

which commutes with the su(l, 1) generators. The superalgebra osp{2, 2) might be constructed by extending su(l, 1) 
algebra with the fermionic generators 

V + = <r + a 2 , V- = o+ai, W + = cr-a+, W- = <J-a 2 . (11) 

The total number operator of the system is given by 

J=~(N-a ). (12) 

The conserved quantity of a physical system possesses osp{2, 2) superalgebra can be written as 

K = N — ~a (13) 

which commutes with the operators of the osp(2, 2) superalgebra. The generators of the osp(2, 2) superalgebra satisfy 
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the following commutation and anti-commutation relations: 

[J+,J-] = -2Jo, [Jo,J±]=±J±, 
[J, J±] = 0, [J, Jo] = 

[Jo,V ± ] = ±\v±, [J ,W±] = ±^W±, 
[J±,V T ] = V ± , [J±,W T ] =w±, 

[J,W±] = -\w±, [J,V±] = \v ± (14) 

[J±,v ± ] = o, [J±,w ± ] = o 

{V±,W ± } = J±, {V ± ,W T } = ±J - J 
{V±,V ± } = {V±,1/ T } = 
{^±,T^±} = {IF±,I^ T } = 0. 

The Hamiltonian of a physical system, with an underlying osp(2, 2) symmetry, has been expressed in terms of the 
operators of the corresponding algebra. In general, the Hamiltonian is exactly solved and the spectrum of the physical 
system can be calculated in a closed form when the Hamiltonian of the system can be written in terms of number 
operator and/or diagonal operator Jo, or it can be diagonalized within the representation [J]. The abstract boson 
and/or fermion algebra can be associated with the exactly solvable Schrodinger equations by using the differential 
operator realizations of boson operators. This connection opens the way to an algebraic treatment of a large class of 
potentials of practical interest [28|. The combinations of the operators of osp(2,2) superalgebra have direct physical 
meaning, being related to the quantum spin systems. 

The Hamiltonian JJJ can be expressed in terms of operators of osp(2, 2) algebra 



H = 2fajJ Q + - k [V+ - V- + W- - W+] + gfiB (N - K) . 



(15) 



where k = y ^-^j^-Xr- Consequently we have shown that the Rashba Hamiltonian possesses the osp(2, 2) symmetry. 

The Hamiltonian (|15|l has (27V + 1) linearly independent eigenfunctions, and it is quasi-exactly solvable. 

Now we turn our attention to the solution of the Hamiltonian (JJJ. A simple connection between the Hilbert 
space and the Bargmann-Fock space can be obtained by transforming the differential realizations of the creation and 
annihilation operators 10. This can be done by introducing the following similarity transformation operators 



A = cxp 



-— (a x a 2 + a 2 ai) 



r = exp 



- (a? + a+ 2 - a\ - 



(16) 



The similarity transformation and change of the variable y — > iy give the following realizations 



h = rAaiA r = 



b 2 = rAa 2 A _1 r _1 = 



m*uj dx ' 
h d 



b+ = TAa+A-^- 1 



6+ = rAaJ A _1 r _1 



-y- 



m*uj dy z 1 V h 

It is obvious that the operator K commutes with the whole Hamiltonian and the eigenvalue problem 



K \nin 2 ) 



k + 2 ) l n i, n 2) ■ 



in the Bargmann-Fock space leads to the following solution: 

iftx, y) = x k (j> (xy) |T) + x k+1 ^{xy) ||) . 



(17) 



(18) 



(19) 



where |f) stands for up state and |J.) stands for down state. The eigenfunction of the Hamiltonian can be obtained 
from the relation 



\ni,n 2 ) = A : r 1 ip(x,y). 



(20) 
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Substitution of (|19fl into the Hamiltonian Q leads to the following set of one variable differential equations 



2z— + k + l + — - - — 
dz 2u> 2hu> hw 



<j>i (z) + X R 



k + 1 



m*ujz d 
h dz 



„ d , „ (k + X)u! c uqB E . 
2z — + k + 2 + V ' + ^— - — U 2 (z + X R 
dz Zlu Znio nw 



m*Lo d 
h dz 



02 0) = o 
0i (z) = 



(21a) 
(21b) 



where z = xy and E is the eigenvalues of the Hamiltonian H and 0i (z) and 2 (z) correspond up and down eigenstates 



i 0i (z) anc 
5is|19j whic 



of the Hamiltonian H, respectively. Following the Reik's analysis 19] which was constructed to obtain the solution of 
the E x e Jahn- Teller Hamiltonian, one can obtain the isolated exact solution of the differential equations. Here we 
present solution of the problem in the framework of the quasi-exactly solvable problem. 



SOLUTION OF THE HAMILTONIAN 

In the previous section we have formulated the Rashba Hamiltonian based on the two boson operators and we have 
discussed its transformation to the one variable differential equation. In this section we present a transformation 
procedure which leads to the quasi-exact solution of the Hamiltonian . The Hamiltonian has been characterized 
by two boson operator. It is possible to transform the Hamiltonian that can be characterized by one boson operator. 
The connection between two and single boson Hamiltonian is given by a similarity transformation induced by the 
metric El 01: 



S = (a+) 



(d_(7 + +(7 + CT_a+) (a+) Q i ai 



The operator (a^") a ^ Ql acts on the state \ni,n 2 ) as follows, 



(22) 



(a+)^ ai |n 1 ,n 2 ) = (a+)" 1 \n u n 2 ) = 



n 2 \ 



(n 2 + ni)\ 

then we can easily obtain the action of the operator S on the two component state 

S(|n 1 ,n a >|T) + |n 1> n a )U» = 
/ (n 2 + ni)\ 



|n 1; n 2 + ni) 



n 2 ! 



f V (ri2 + ni + 1) \ni,n 2 + n\ + 1} |f) + \ni,n 2 + n{) |J.)^ 



(23) 



Since a\ and a 2 commute, the transformation of a\ and under S can be obtained by writing 
[ ai ,b} = [at,b] = 0, 



(24) 



with 



S 1 aiS — aia^; S 1 afS — [a 2 ) 1 



and transformation of a 2 and at as follows 



(25) 



S~ 1 a 2 S = a 2 + (a^ai + er + er_)(aj) _1 ; S ,_1 a 2 f S = ■ 
Similarly the transformation of the Pauli matrices are 

S~ 1 a + S = <J+{a2 S~ 1 <T-S = a-ci^; S~ 1 aoS = <jq 

The Hamiltonian H can be transformed 

H = S^ 1 HS = hui(2a^ai + a\a 2 + (7 + <r_ + 1) ^~ {a 2 a 2 + c_|_cr_ 

k [(1 — a 1 )a + + (a 2 a 2 f + afai + <7 + er_ — a+ )er_] + -g[iBa 



(26) 



(27) 



(28) 
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It is obvious that H can be characterized by a fixed number a\ai = — j — 1. Here j takes integer values. The 
transformed Hamiltonian H includes one boson operator and possessing infinitely many finite-dimensional invariant 
subspaces with a basis function: 

-(£$))• (29) 

where 4>{z) is two component spinor and p n (z) and q n +\(z) are polynomials of degree n and n + 1 respectively. 

a;, we 



Substituting Bargmann-Fock space realizations of the bosonic operators l|17l) by changing variable z = y —j— 
obtain the following single variable differential equation: 



H 



2z- j + (T + cr_ 

dz 



hljj r 



CO 



f 



d 

dz 



(7+ + cr_ (z 



dz 



The eigenvalue problem can be expressed as 

H(/)(z) = E<f>(z) 

The action of the H on the basis function <fi(z) gives the following recurrence relation: 

[2nfvjJ-E + e j +e b }p n {E)-K[q n+ i{E)-{n+l)q n {E)] = 
[2nhuj - E + - e b ]q n+1 (E) + K [p n+1 (E) + {j - n)p n (E)} = 



(30) 



(31) 



(32) 



where p n {E) and q n (E) are the coefficients of the polynomials p n (z) and q n (z), respectively. The justified parameters 
are given by 



fkt> 



1 . (3 + 2j> t 

2" J + ^T- 



(33) 



It is necessary that the determinant of these sets of be equal to zero giving the compatibility conditions that establish 
the locations of the exact eigenvalues on the energy baseline. The recurrence relation implies that the wavefunction is 
itself the generating function of the energy polynomials. If E% is the roots of the the eigenfunction truncated at 
n = j and it is the exact eigenvalues of the Hamiltonian H. These recurrence relations can be written in the matrix 
form: 



E+ 


— K 















Po 


JK 


E_ 




K 










H 





2 k 


E+ 


+ 2hw 


— K 







Pi 








(j 


- 1)k E„ 


+ 2hu) 


K 




12 













3 k 


E + +3huj ■ 




P2 



= 



(34) 



where E± = ej ± — E. It is obvious that the determinant of the matrix forms polynomials in E and the first few of 
them are given by 



D {E) = E+E_ 

Di(E) = (E- + 2hw) [E+E-(E+ + 2hw) - n 2 (E+ - 2hw)} 
D 2 (E) = (E- + Ahw) [E+E-{E+ + 2hu>)(E_ + 2fku)(E + + Ahw) 

+ n 2 (2hwE + (E- + Ahw) + lQ2h 2 uj 2 {E- + 2hw) - 2E\EJ) 

+ 2k 4 (E + - 2hw)} 



(35) 



for the values j = 0, 1 and 2, respectively. The eigenvalues E of the Hamiltonian Q can be determined by finding 
the roots of the polynomials (|35|) . 
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CONCLUSION 



We presented a QES to the problem of an electron in a quantum dot in the presence of both the magnetic field 
and spin-orbit coupling. Our formulation gives a unified treatment of the solution of the spin-boson physical systems. 
It has been shown that the mathematical structures of the quantum dot Hamiltonian including Rashba spin-orbit 
interaction and some of the spin-boson physical systems are identical. We have also shown that the Hamiltonian 
possesses osp(2, 2) hidden symmetries. The suggested approach can be modified to solve the quantum dot Hamiltonian 
including Dresselhaus interaction. 

Furthermore we have presented a transformation procedure that offers several advantageous, especially if one wishes 
to describe the eigenvalues of the Hamiltonian. It is obvious that the algebraic techniques have been used in a variety of 
problems to compute their spectrums. We have presented the steps towards an extension of the algebraic formulation 
of the quantum dot Hamiltonians. 

The technique given in this article can be extended in several ways. The Hamiltonian of a quantum dot including 
position dependent effective mass may be formulated and solved within the procedure given here. We hope that our 
method leads to interesting results on the spin-orbit effects in quantum dots in future research. 



* Electronic address: tutunculer@gantep.edu.tr 

* Electronic address: koc@gantep.edu.tr 

* Electronic address: olgar@gantep.edu.tr 

[1] Bychkov Y A and Rashba E I, 1984 J. Phys.C 17 6039 
Datta S and Das B 1990 Appl. Phys. Lett. 56 665 
Winkler R 2000 Phys. Rev. B 62 4254 
Prinz G A 1998 Science 282 1660 

Wolf S A, Awschalom D D, Buhrman R A, Daughtan J M, vov Molnar S, Roukes M L, Chtchelkanova A Y and Treger D 
M 2001 Science 294 1488 
Dresselhaus G 1955 Phys. Rev. 100 580 

Wang X F, Vasilopoulos P, Peeters F M 2002 Appl. Phys. Lett. 80 1400 
Bandyopadhyay S 2000 Phys. Rev. B 61 13813 
Governale M 2002 Phys. Rev. Lett. 89 206802 

Valin-Rodriguez M, Puente A and Serra L 2003 Eur. Phys. J. B34, 359 

Tsitsishvili E, Lozano G and Go golin A O 2003|htt'p://arxiv.org/cond-mat/0310024| 
Frustaglia D and Richter K 2003 http://arxiv.org/cond-mat/0309228 
de Sousa R and Sarma D S 2003 Phys. Rev. B 68 155330 
Turbiner A V and Ushveridze A G 1987 Phys. Lett. A 126 181. 
Bender C M and Dunne G V,1996 J. Math. Phys. 37, 6 

Gonzales- Lopez A, Kamran N and Olver P J 1993 Comm. Math. Phys. 153 117. 
Shifman M A 1989 Int. J. Mod. Phys. A4 2897. 
Judd B R 1979 J. Phys. C 12 1685. 

Reik H G, Stiilze M E and Doucha M 1987 J. Phys. A 20, 6327. 
Loorits V 1983 J. Phys. C 16, L711. 
Pooler D R 1978 J. Phys. A: Math. Gen. 11 1045 

Kog R, Tiitiinriiler H, Koca M and Korciik E 2003 Prog. Theor. Phys. 110 399 
Jaynes E T and Cummings F W 1963 Proc. IEEE 51 89 
TurEA 2000 Opt. Spectrosc. 89 574 

Jing-Bo Z U and Xu-Bo Z O U 2001 Chin. Phys. Lett. 18 51 

Chen Yong-Qing 2000 J. Phys. A: Math. Gen. 33 8071; 2001 Int. J. Theor. Phys. 40 1113; 2000 Int. J. Theor. Phys. 39 
2523 

Chen Yong-Qing, Xiao-Hui Liu and Xing-Chang Song 1994 Commun. Theor. Phys. 22 123 
Alhassid Y, Giirsey F and Iachello F 1983 Ann. Phys. (N. Y.) 148 346 
Tiitiinciiler H and Kog R 2004 Pramana J. Phys. 62 993 



